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Abstract 

The inflow problem of full compressible Navier-Stokes equations is considered 
on the half line (0,+oo). Firstly, we give the existence (or non-existence) of the 
boundary layer solution to the inflow problem when the right end state (p+, u+, 9+) 
belongs to the subsonic, transonic and supersonic regions respectively. Then the 
asymptotic stability of not only the single contact wave but also the superposition 
of the boundary layer solution, the contact wave and the rarefaction wave to the 
inflow problem are investigated under some smallness conditions. Note that the 
amplitude of the rarefaction wave can be arbitrarily large. The proofs are given by 
the elementary energy method. 

1 Introduction 

In this paper, we consider the half space problem of the full (or non-isentropic) compress- 
ible Navie-Stokes equations in Eulerian coordinate: 



Pt + (pu).^ = 0, X > 0, t > 0, 

{pu)t + {pv^ + p)x = pu^x, X > 0, t > 0, 

[p(e + it)] + [pu{e + — ) + = i^S^^ + p{uUx)x, x > 0, t > 0, 
Z t I 



where p{t,x) > is the density, u{t,x) is the velocity, 9(t,x) is the absolute temperature 
of the gas, and p = p{p, 6) is the pressure, e = e(p, 0) is the internal energy, /x > is the 
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viscosity constant, and k > is the coefficient of heat conduction. Here we consider the 
perfect gas, that is 



p = Rp9 = Ap'^ exp ( — =— s), e = \-const., (1.2) 

it 7 — 1 

where s is the entropy, 7 > 1 is the adiabatic exponent, and A,R>0 are gas constants. 
The initial values are given by 

{p,u,9){t = 0,x) = {po,uo,Oo){x) {p+,u+,9+),x ^ +00, (1.3) 

where (p+, m+, 6'_|_) is a constant state with p+,0+ positive. The boundary values are the 
following: 

(p,«,^)(t,a: = 0) = (p_,M_,^„), (1.4) 

where p_ > 0, 6'„ > 0, m_ are given constants. And of course the initial values (11. 3p and 
the boundary condition (11.41) satisfy the compatible condition at the origin (0,0). 

According to the sign of the velocity m_ on the boundary {x = 0}, the following three 
types of problems are proposed [T7]: 

(1) the inflow problem, i.e., the velocity U- > 0; 

(2) the outflow problem with < 0; 

(3) the impermeable wall problem, i.e., U- = 0. 

It should be remarked that in the cases (2) and (3), the density p_ can not be given on 
the boundary by the theory of well-posedness on the hyperbolic equation (ll.ll) i . 

In this paper, we are interested in the case of the inflow problem (II. ip . (I1.3p - (I1.4I) . 
When K = p = 0, the compressible system (II. ID becomes the inviscid Euler system 

Pt + ipu)x = 0, 
{pu)t + {pu^ +p)x = 0, 

W + v)] + [^"(^ + + 

It z 

The Euler system (11.51) is a typical example of the hyperbolic conservation laws. It is 
well-known that the main feature of the solutions to the hyperbolic conservation laws is 
the formation of the shock wave no matter how smooth the initial values are. The Euler 
system (II. 5p contains three basic wave patterns in the solutions to the Riemann problem. 
They are two nonlinear waves, called shock wave and rarefaction wave, and one linear 
wave called contact discontinuity. The above three dilation invariant wave solutions and 
their linear superpositions in the increasing order of characteristic speed, i.e., Riemann 
solutions, govern both local and large-time behavior of solutions to the Euler system. The 
invscid Euler system (II. 5p is an ideal model in gas dynamics when the dissipation effects 
are neglected, thus it is of great importance to study the corresponding viscous system 



(1.5) 

= 0. 
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There has been a large hterature on the large-time behavior of the solutions to Cauchy 
problem of the compressible Navier- Stokes equations (11.11) toward the viscous versions of 
the three basic wave patterns. In 1985, Matsumura-Nishihara [H] firstly proved the 
stability of the viscous shock wave to the isentropic compressible Navier-Stokes equations 
(i.e., the entropy s is assumed to be constant and the energy conservation law is not 
considered). Since then, many authors had been attracted to study the stability of the 
viscous wave patterns and much progress has been made. We refer to [6], [7], [9], [11], 
[T2] . [13], [H], [15], [19], [21], [25], [26] and some references therein. All these results show 
that the large-time behavior of the solutions to Cauchy problem are basically governed 
by the Riemann problem of the corresponding Euler equations. 

Recently, the initial-boundary value problem (IBVP) of (11.11) attracts increasing inter- 
est because it has more physical meanings and of course produces some new mathematical 
difficulties due to the boundary effect. Not only basic wave patterns but also a new wave, 
which is called boundary layer solution (BL-solution for brevity) [T7], may appear in the 
IVBP case. Matsumura [17] proposes a criterion on the question when the BL-solution 
forms to the isentropic Navier-Stokes equations. The argument is also true to the full 
Navier-Stokes equations (11.11) . Consider the Riemann problem to the Euler equations 
(II. 5p . where the initial right end state is given by the far field state (p+, 9+) in (11.31) . 
and the left end state (p_, u^,9J) is given by the all possible states which are consistent 
with the boundary condition (11.41) at {x = 0}. Note that for the outfiow problem, p_ 
can not be prescribed and is free on the boundary. On the one hand, when the left end 
state is uniquely determined so that the value at the boundary {x = 0} of the solution 
to the Riemann problem is consistent with the boundary condition, we expect no BL- 
solution occurs. On the other hand, if the value of the Riemann problem's solution on 
the boundary is not consistent with the boundary condition for any admissible left end 
state, we expect a BL-solution which compensates the gap comes up. Such BL-solution 
could be constructed by the stationary solution to Navier-Stokes equations. The existence 
and stability of the BL-solution (to the infiow or outfiow problems, to the isentropic or 
full Navier-Stokes equations) are studied extensively by many authors, see [2], [1], [H], 
[16], [Uj ^20j, [22j, [27j, etc.. For the infiow problem of the full Naier-Stokes equation 
(ll.ip - (ll.4l) . Huang-Li-Shi [2] proved the stability of the BL-solution in some cases. More 
precisely, they show that when {p±,u±,6±) both belong to the subsonic region, the BL- 
solution is expected and the stability of this BL-solution and its superposition with the 
3-rarefaction wave is proved under some smallness assumptions. Notice that both the 
BL-solution and the rarefaction wave must be weak enough. When the boundary value 
(p_, M_, 9-) belongs to the supersonic region, there is no BL-solution. Thus the large-time 
behavior of the solution is expected to be same as that of the Cauchy problem and the 
stability of the rarefaction waves is given. 

In this paper, firstly we give the existence (or non-existence) of the BL-solution when 
the right end state (p+, m+, 9+) belongs to the subsonic, transonic and supersonic regions, 
respectively. The rigorous proof is given in Appendix. Notice that it is more natural to 
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present the classifications according to the locations of the right end state (p+, m+, 6'+) 
from the qualitative theory of the autonomous ODE system. Then we prove the stability 
of not only the single contact wave but also the superposition of the BL-solution (subsonic 
case), the viscous contact wave and the 3-rarefaction wave to the inflow problem (11 .ip - 
(ll.4p . Here the amplitude of the rarefaction wave can be arbitrarily large. 

Now we briefly review some key analytic techniques in studying the stability of the 
basic wave patterns. The strict monotonicity of the corresponding characteristic speed 
along the wave profiles plays a crucial role in stability analysis of the viscous shock wave 
and rarefaction wave. Precisely speaking, the shock wave is a compression wave, so the 
characteristic speed is monotone decreasing in the shock profile. Thus anti-derivative 
variable to the perturbation should be introduced in the stability analysis. While the 
rarefaction wave is an expansion wave and the characteristic speed is monotone increasing 
along the rarefaction wave, thus the direct energy estimates to the perturbation itself are 
available. However, the characteristic speed along the contact wave is constant, and the 
spatial derivative of the velocity changes signs along the contact wave profile. Due to the 
degenerate characteristics, the stability of the contact wave profile to the compressible 
Navier-Stokes system (ll.ip is just proved by |6j and j9j in 2005, twenty years later than 
the nonlinear wave in 1985. In [6] and [9], the anti-derivative variable to the perturbation 
is introduced and the proof framework is motivated by the viscous shock profile. Notice 
that a convergence rate of the order of (1 + t)~4 in sup-norm is a by-product of the 
estimation. However, there is no convergence rate obtained so far for the viscous shock 
wave and the rarefaction wave. 

Recently, Huang-Matsumura-Xin [7J obtained a new estimate on the heat kernel which 
can be applied to the study of the stability of the viscous contact wave in the framework 
of the rarefaction wave, see [3] or Lemma 3.4 in the present paper. Namely, the anti- 
derivative variable of the perturbation is not needed and the estimations to the pertur- 
bation itself are also suitable to get the stability of the viscous contact wave. But the 
time-decay rate can not be gotten as a compensation. More importantly, the advantage 
of this framework is that it can be used to study the stability of the contact wave to the 
IBVP of (II. ip since the boundary terms could be treated conveniently. We will make full 
use of this new estimate on heat kernel to study the inflow problem (II. ip and get the 
expected results. 

The novelty of the paper lies in the following three aspects: (1) The rigorous proof 
and the classifications of the existence (or non-existence) of the BL-solution to the inflow 
problem. (2) The stability of the superposition of three different wave patterns (the BL- 
solution, the viscous contact wave and the rarefaction wave). (3) The large amplitude of 
the rarefaction wave in the superposition wave. The main difficulties in our proofs are 
how to deal with the boundary terms and the interactions of three different wave patterns. 

Because the system (II. ip we consider is in one dimension of the space variable x, it is 
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convenient to use the following Lagrangian coordinate transformation: 

p{y,t)dy, t^t. 



X 



Thus the system (11. ip can be transformed into the following moving boundary problem 
of Navier-Stokes equations in the Lagrangian coordinates: 



vt-u^ = 0, 



Ut+Px= K — )x, 

(e + —)t + [pu)x = n{—)x + /i( )x, 

2 V V 



;i.6) 



X > cr_t, t > 0, 
X > cr^t, t > 0, 

X > a^t, t > 0, 

1 

(f , M, 6')(t, X = (T_t) = (f_, 6*-), ti_ > 0, 

, {v,u,9){t = 0,x) = {vo,uo,9o){x) ^ {v+,u+,9+), as x ^ +oo, 

where v(t,x) = represents the specific volume of the gas, and the boundary moves 
with the constant speed cr_ = — — < 0. Now we have that for the perfect gas, 



V R 



R 



7-1 



-6 + const. 
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In order to fix the moving boundary x = cr^t, we introduce a new variable ^ = x — a^t. 
Then we have the half space problem 



Vt — cr_fg — = 0, 



V 

ie+^ t-^- e+— 5 + Mc = '« - 5 + /" — S, e>0, t>0, 
2 2 V V 



e>o, t>o, 
e>o, t>o, 



■u_ > 0, 

as ^ — s> +00. 



Given the right end state we can define the following wave curves in the 

phase space (f , u, 9) with v > and 9 > 0. 
• Contact wave curve: 



CD{v+,u+,9+) = {{v,u,9)\u = u+,p = p^,v^ v+}. 
BL-solution curve (subsonic case, i.e., m+, 6'+) G ^l^ub)'- 



;i.9) 



BL{v+,u+,9+) 



-a_ = ^,{u,9)eM{u+,9+)], (1.10) 



where A^(m+, 9+) is the center-stable manifold defined in Lemma 2.1 below. 
• 3-Rarefaction wave curve: 



i?3(t;+,M+,6'^ 



V > v+, u = — 



X-i{r],s+)dr], s{v,9) 



1.111 



where 5+ = s(f+, 6'+) and A3 = A3(f , s) is the third characteristic speed given in ( 12. ip . 
Our main stability results are, roughly speaking, as follows: 

(I) . If the state m_, 9-) e CD(f+, u+,9+), then the viscous contact wave is asymp- 
totic stable under some smallness conditions which are given in Theorem 2.1. 

(II) . If the state (f _, 6'_) G BL-CD-R3(w+, u+, 6*+), then there exist a unique state 

9*), such that (f_,M_,6'_) G BL(f*, m*, 6**), 



6**) G and a unique state 



{v^,u^,9^) G CD{v*,u*,9*) and (t;*,M*,6'*) G R3(f+, m+, 61+) and the superposition of the 
BL-solution, the viscous contact wave and the rarefaction wave is asymptotically stable 
provided that |(m_ — u^,9_ — 9^:)\ and |f* — v*\ are suitably small and the conditions in 
Theorem 2.2 hold. It is remarked that the BL-solution and the viscous contact wave must 
be weak but the rarefaction wave is not necessarily weak. 

Notations: Throughout the paper several positive generic constants are denoted by 
c, Co, Ci, ■ ■ ■ or C, Ci, C2, ■ ■ ■ without confusions. The small constant u > is used in Young 
inequahty 

ab < vqP^ + CJF\ — + — = 1, 



1 1 



where is the constant depending on v. For functional space, i7'(R+) 
/-order Sobolev space with the norm 



denotes the 



where||/|| 



2 Preliminaries and main results 

In this section, we will show the wave patterns considered in the paper. We start with 
the BL-solution. 



2.1 BL-solution 

It is known that the hyperbolic system (II. 5p has three characteristic speeds 

A, = 





(2.1) 
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The sound speed c{v, 0) and the Mach number M are defined by 



c(v,e) =vj— = Jr^, (2.2) 
V V 

and 



Miv,^,0) = ^\, = ^=, (2.3) 
c{v, 9) ^/R'y6 

respectively. 

We divide the phase space {{v,u,9),v > 0,9 > 0} into three regions: 



sub 



n 
p 

trans 
^ super 



{iv,u,9)\Miv,u,9)<l}, 

{{v,u,9)\M{v,u,9) = l}, (2.4) 
{{v,u,9)\ M{v,u,9) > 1 }. 



Call them the subsonic, transonic and supersonic regions, respectively. If adding the 

^sub^ 



alternative condition m > or m < 0, then we have six connected subsets fi^^, Tf^ans^ ^ind 



super' 

When (f_,M_,6'_) G ^~^uby we have 

Ai('U_,M_,6'_) < (T- < 0, 

hence the existence of the traveling wave solution 

{V^,U^,e^m = {v.,u.,9.), (\/^f/^,e^)(+oo) = (t;+,«+,e+), 

to (11. 6p . or the stationary solution to (II. 8p is expected. We call this traveling wave solution 
{V^ , U^, Q^){C,) the boundary layer solution to the infiow problem (11.61) . Note that the 
speed of the traveling wave is just the speed of the moving boundary of (II. 6p . 

In the following, we will give the existence (or non-existence) of BL-solution to the 
infiow problem (II. 6p . From (12.51) . BL-solution [V^ satisfies the following 
ODE system 

-a^v^y - {u^y = 0, e>o, 
-a_(t/^)' + (p^)' = M^)', e>o, 

I (V^,f/^,e^)(0) = (t;_,w_,^^_), {V^U^,Q^){+^) = {v^,u+,9^). 

(2.6) 

where :=p(y^,e^) - ^^"^ 



Integrating the system (12.61) over +00) implies that 



/i 



[ (t/^e^)(o)=/«_,^.), (f/^,e^)(+oo) = K,^^+). 



Let ^ = in fl2.7P i . we have 



which is the first condition of BL-solution curve in (11 .101) . 

From the fact m_ > and f± > 0, we find that m+ must satisfy 



u+ = — w+ > 0. 



Rewrite (I22D2-(EZD4 as 



9, 



-V' 



fi:(7 — Ij K 2k 



Denote 



and 



J 



ul-Re+ R 

Re+ Ru+ 
K k(7— 1) 



f (M^7-l)u+ ^ 



where M+ = M(t;+, m+, 6*+). 

Then we obtain the automatons ODE system 



QB 
tB t^B 



J 



e 



B 



+ 



(f/^, e^)(o) = (n_ - - (f/^, e^)(+oo) = (0, 0) 



where 



Re, 



B\2 



U 



7^ (U 



B\2 



+ 



R 



- 1) 



(2.14) 



Now we state the existence results of the solution to fl2.13p while its proof will be shown 
in Appendix. 

Lemma 2.1 (Existence of BL-solution) Suppose that v± > 0, U- > 0, 6± > and 

let 6^ = |(m+ — n_, 6^+ — 9^)\. If < 0, then there is no solution to ^2.1Ch) or Ii2.13\) . If 
> 0, then there exists a suitably small constant 6 > such that if < 6^ < 6 , then 
Case I. Supersonic case: M+ > 1. Then there is no solution to (2.10) or (2.13). 
Case II. Transonic case: M+ = 1. Then {u^,9^) is a saddle-knot point to Ii2.13\) . 

Precisely, there exists a unique trajectory T tangent to the line 



/iM+(f/^-M+)-K(7-l)(e 



B 



9, 







at the point (m_|_, 9^ 
satisfying 



B 



For each {u_,9_) G T, there exists a unique solution {U^,Q^) 



< c- 



n 



0,1,2, 



(2.15) 



Case III. Subsonic case: M+ < 1. Then the equilibrium point ('U_(.,^+) is a saddle 
point of Ii2.13\) . Precisely, there exists a center-stable manifold M. tangent to the line 

(1 + a2C2U+){U^ - u+) - 02(6^ 



on the opposite directions at the point (ii+, 9 
equation 



MlR-f k(7 - 1) 



Here C2 is one of the solutions to the 



y 



MlRjK 







and a2 



R 



with > 0, < are the two eigenvalues of the matrix A. Only 



when {u-,9^) G A4, does there exist a unique solution {U^,Q^) C M. satisfying 



Jn 



n 



0,1,2, 



(2.16) 



Remark: This Lemma is the first one for the classifications of the BL-solution to the 
inflow problem (1.1). The stability of the single BL-solution in subsonic case (Case III) 
is proved in [2]. In this paper, we are concerned with the stability of the superposition 
of this subsonic BL-solution with viscous contact wave and rarefaction wave. As for the 
stability of the BL-solution in transonic case (Case II) and its superposition with other 
wave patterns, it is in consideration [23j . 
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2.2 Viscous contact wave 

If {v_,u_,9_) e CD{v+,u+,9+), i.e., 

M_ = u+, p„ = p+, (2.17) 
then the following Riemann problem of the Euler system 

Vt-U:r: = 0, 

Ut+p^ = 0, (2.18) 

with Riemann initial value 

(f , u, 9)(t = 0,x 
admits a single contact discontinuity solution 

{v,u,e){t,x) = 



(f M_, 6'_), X < 



(■u_,M_,6'_), X < 0, t > 0, 
(t;+,M+,6'+), X > 0, t > 0. 

From [6] , the viscous version of the above contact discontinuity, called viscous contact 
wave [V^^, If^^, Q'^^){t,x), could be defined by 

e^^(t,x) = e^*'"( ^ ^ 



p+ 



V''''{t,x) = (2.19) 

U^^itx)-u I ^(7-i)eg^(t,x) 

where B'^*'^(?7), rj = — ^==, is the unique self-similar solution of the nonlinear diffusion 
equation 

_ Kp+il - 1) f 
ft — 



i?27 V©/.' (2-20) 

e(±oo) = 

Thus the viscous contact wave defined in f l2.19p satisfies the following property 

|e^^ - ^±1 + (1 + t)i|e,^^| + (1 + t)|e^i^| + (i + t)i |e,Tj = o(i)5^^e- (2.21) 

as |x| — >■ +00, where 6'"^ = — (^-\ is the amplitude of the viscous contact wave 
and Co is a positive constant. Note that ^ = x — (T_t, then (V'-^^ ,U^^ ,Q^^){t,x) = 



10 



^yCD^ jjCD^ qCd^^^^ ^ ^ Satisfies the system 



tCD 

R 



Vc 



c 



e 



CD 



yCD 



+ Q2 



1 T~>cn RQ'^^ 
where P = -yen' = 



p_ and the error terms Qi, Q2 are given by 



3 CQ(i + a_t)^ 



as 1^ + a_t| — » +00, 



Q2 



s:CD\2 



o(i)(|ef^r + |effp) 



0(l)(5^")'(l + t)-'e 



eoa + .T_t)^ 
1 + t 



as 1^ + cr_t| — s> +00. 



(2.22) 



(2.23) 



2.3 Rarefaction wave 

If {v-,U-,9-) G Rs{v+,U-^,9+), then there exists a 3-rarefaction wave {v^ , u"^ , s^) {x / 1) 
which is the global (in time) weak solution of the following Riemann problem 



r r r\ 

v^-u^ = 0, 



t > 0,x e R, 



7-1 
(^;^n^e^)(0,x) 



(2.24) 



(f_, OJ),x < 0, 

(f+,M+,^+),X > 0. 



From [5], it is convenient to construct the approximated rarefaction wave (V^, U^, G^)(t, x) 
to the infiow problem (11.61) by the solution of the Burgers equation 



wt + wwx = 0, 
w{0, x) = wo{x) 



w_, X < 0, 

/•ex 

w^ + Ca6' / y'^e-ydy,x > 0, 



(2.25) 



where 



g > 16 is some fixed constant, Cg is a constant such that 



Cq y'^e ^dy = 1, and e < 1 is a small positive constant to be determined later. The 
solution of the Burgers equation (12.251) have the following properties: 
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Lemma 2.2 ([5]) Let < < Burgers equation i\2.25^ has a unique smooth 
solution w(t,x) satisfying 

i) W- < w{t,x) < Wxit,x) > 0, 

ii) For any p {I < p < oo), there exists a constant Cpg such that 

II w,{t) |Up< Cp,mm{6re'-'/P, SI/H-'^'/^}, 
\\w,,{t) \\L,<Cp,mm{6re^-'/P, {61/^ + 6l/'^)t-'-^'/'^}, 

hi) If X < W-t, then w{t, x) = W-, 

iv) sup\w(t,x) — w^{x/t)\ ^ 0, as t ^ oo. 

Thus we construct the approximated rarefaction wave {V^, U^, Q^){t, x) by 



5^ = s(V^,e^) = s+, 

w{l + t,x) = XsiV''{t,x),s+), 



U^{t,x) = u+- j 

J V-\. 



(2.26) 



X'iiv, s+)dv. 



Note that ^ = x — cr_t, then the smoothed 3-rarefaction wave iV^, U^, G^)(t, ^) defined 
above satisfies 



Vt"" - a.Vf - f/« = 

tR „ ttR 



- a^Uf + P^ = 0, e > 0, t > 0, 



^(ef - a_ef ) + p^?7f = 0, ^^-^^^ 

[ (V,f/^e^)(t,O) = (t;_,«_,0_), (r^[/^e^)(t,+oo) = (t;+,«+,^+), 
where := p(\/^, 6-^). 

Lemma 2.3 ([5J) Let 6^ = |(f+,M+,^+) — (w_,m_,^_)|. The approximated rarefaction 
wave {V^, f/^, e^)(t, satisfies 

i) Fore>0, ,t>0, f/f(t,O>0, 

ii) For any p (1 < p < oo), there exists a constant Cpg such that for t > 0, 

II (l^f,[/f,ef)(t) |Up< Cp,min{5V-i/^ (5^)^/P(l + t)-i+i/n, 

II iV^lU^%Qf^m |Up<Cp,min{5V-i/P, ((5^)1/^ + (5«)i/^)(l + t)-i+i/n, 

iii) If ^ + (T_t < A3(t;_,M_,^_)(l + t), then (\/^, ?7^, e^)(t, = {v-,u_,e_), 

iv) sup I {V^, f/^, e^) (t, - (^^^ m", ^0 (ifc) I ^ 0' as t ^ oo. 



2.4 Main results 

Now we can state our main results. The first one is in the following. 
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Theorem 2.1 If G CD{v+, u+, 6+). Let {V^^ , U^^, e^^)(t, x) be the viscous 

contact wave defined in (12.191) . There exists a smaU constant 6^ such that if the wave 
amphtude 5*"^ and the initial values satisfy 

then the moving boundary problem (II. 6p or the half space problem (II. Sp admits a unique 
global solution {v,u,9){t,^) satisfying 

(t;- V^^,M-f/^^,^^-0^^)(t,O G C([0,+oo),i/i(R+)), 
(^_yCD)^(t,^)eL2(0,+oo;L2(R+)), 



\CD 



and 



lim sup |(t;-\/^^,u-[/^^,^-0^^)(t,OI =0. 



Now we state our second result. If G BL-CD-R3(f+, u+, then there 

exist states G O^j, and (f*,M*,6'*), such that (i;„,m_,6'_) G BL(t>*, m^,, 6'*), 

(f^,u^,6'*) G CD(f*, n*, 6^*) and {v*,u*,9*) G R3(i;+, m+, 61+). In fact, by(f _, m_, 6'_) G 
BL(t>^,, M^, 6'*) and (ll.lOp . we have 



and by (w*,^*,^) G R3(i;+, u+, ^+), ([m]) gives 



\3{r],s+)dr], V* > v+. 



(2.28) 



(2.29) 



Thus the two curves (I2.28P and (12.290 have a unique intersection point m = -u in {v, u) 
space. If = u* = u, then = v*, thus there is no contact wave. By {v^,u^:,9^:) G 
CD{v*,u*,9*), we have 

Thus if = u* ^ u, then v^, ^ v* and there exists a contact wave. Among the three 
values M*(= M*), and f*, only one is independent, the other two can be determined 
accordingly. 

Now assume that u*(= u*) is given, then from (I2.28p and (12.290 . we can determine 
and f * by 



V = Va 



7-1 



:(-U* - «+) + 1 



2 

1-7 



2A^R-i9^ 

By the definition of the rarefaction wave curve -R3 in (11.110 . we have s{v^,, 9^) = s_|_, i.e., 

2 



9* = 9, 



9, 



7-1 
2Av/i?7^H 



(m* - + 1 
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Again by the contact wave curve (11 .Qp . one must have 



P* = P*, 



thus 



e* 



e. 



7-1 



27 
7-1 




'+ 



(m* - M+) + 1 



(2.30) 



So if large enough, then (n*, 6**) in (12.301) must belong to the region := {(m*, 6'*)|0 < 
< y/Ry9^}. Moreover, from the definition of the BL-solution, (m_,6'_) G J^{u^,9^). 

Thus (u*, 6'*) can be determined uniquely if (m_, 6'_) is given suitably. 
Define the superposition wave (V, U, 9)(t, ^) by 



where {V^ ,Q^){t,^) is the subsonic BL-solution defined in Lemma 2.1 (Case II) 
with the right state (f+,M+,6'+) replaced by {v^,u^,9^), {V^^,U'^^,Q'^^){t,^) is the 
viscous contact wave defined in (12.191) with the states (f_, m_, 9^) and {v+, u^,9^) replaced 
by (v*,M*,6'*) and {v*,u*,9*), respectively, and (y^,U^,Q^)(t,^) is the smoothed 3- 
rarefaction wave defined in (12.261) with the left state (f_,u_,^_) replaced by {v*,u*,9*). 

Theorem 2.2 If {v_,u.,9_) e BL-CD-Rs{v+, u+,9+). Let {V,U,e){t,x) be the super- 
position of the BL-solution, the viscous contact wave and the rarefaction wave defined in 
(I2.3ip . There exists a small constant such that if the BL-solution amplitude 6^, the 
contact discontinuity amplitude 6'-^^ and the initial values satisfy 



then the inflow problem (11.61) or the half space problem (11.81) admits a unique global 
solution {v,u,9){t,^) satisfying 




(2.31) 



'0, 



{v -V,u-U,9 -Q){t,i) e C{[0,+oo),H\R+)) 

{v - V)^{t,0 e L\0,+oo;L\R^)), 

{{u - U)^, {9 - e)^){t, e L^O, +00; H\R+)), 



and 



lim sup \{v-V,u-U,9-Q){t,^)\=0. 
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3 Stability analysis 



In this section we will prove our main stability results Theorem 2.1 and Theorem 2.2. We 
will focus on the proof of Theorem 2.2, i.e., the stability of the superposition wave. The 
proof of Theorem 2.1 is almost same as Theorem 2.2 and we will omit it for brevity. 

Besides the intrinsic properties of the BL-solution, the viscous contact wave and the 
rarefaction wave in the stability analysis, the interaction between the wave patterns should 
be dealt with carefully in the stability analysis. Here we will use the elementary energy 
methods to prove Theorem 2.2 by the classical continuum procedure. 

Firstly we will reformulate the system of the superposition wave (V, U, 9)(t, ^) defined 
in dMH). 



3.1 Reformulation of the problem 



Recall the definition of the superposition wave (V, U, Q){t,^) defined in (12.311) . Then we 
have 



0, 



R 



e>o, t>o. 



u. 



-(6, - a^e^) + PU^ = k{^)^ + /i^ + Q2 



7 



(3.1) 



t (V,f/,e)(t,0) = (t;_ + V^^-t;.,M_ + [/^^-M„^_ + 0^^-^,)(t,O), 
where P = p{V, 0) = and the error terms Qi (i = 1, 2) are given by 



Qi 

Q2 



{PU^ 



pCD _ pH)^ _ ^ 



CDttCD 



P^^'U, 



V 



P'^U?) - K 



yCD>^ 



V 



yCD 



+ Q2. 



• V 



,ef. 



yCD' 



(3.2) 



and Qi 



1, 2) are the error terms defined in (I2.23P to the viscous contact wave. 



Due to the different propagation speeds of the BL-solution, the viscous contact wave 
and the rarefaction wave, we can get the following estimates of the error terms Qi{i = 1, 2): 



R\ 



Q, = o(i)[|(f/f,yf,ef,[/f^)||(y-v^,e-e^,yf^,f/f^,v^f,f/| 

+ \{v^'',Qfmv-v^,Q~ e^)| + |(f/f^,f/f\/f )|] + IQi 

+ \Qi\, 



(3.3) 



R -irR\l2\ 



0(1)(5^ + 5^^)e-'=^i^i+'^ + 0(l)(|f/^^|, |(f/f , 
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for some positive constant c independent of ^ and t. Similarly, 

Q2 = o(i)(5^ + 5^^)e-^«+*) + o(i)(|eg|, |(ef , , f/f )n + IQ2I. 

Denote the perturbation by 

(0, ^, c)(t, = (v, u, 9){t, - {V, u, e)(t, 0, 

then we have the initial boundary value problem of the perturbation {(j),ip,(){tyO'- 



(3.4) 



0t - cr_05 - = 0, 

- a_ij^ + ip- P)^ = - - Qi, 

^(C, - a_Cc) + - PU^) = /.(^ - + /.(^ 
(0, ^, C) (t, e = 0) = - \/, M_ - f/, - 6) (t, e = 0), 
(0, C)(t = 0, = (00, V'o,Co)(0^ (0,0,0), as e^+cx). 



^ > 0, t > 0, 
^ > 0, t>0, 

)-Q2, e>o, t >o. 



(3.5) 

Since the local existence of the solution of (13.51) is well-known, we just state it and 
omit its proof for brevity. 
Denote that 

N{t)= sup ||(0,^,C)(r,-)||?, 



and define the solution space by 



(0,^,C)(t,O 



re[o,t] 



(0,^,C)(t,OeC([O,T];i7i(R+)), 

(^5,C5)GL2(0,T;ifi(R+)), 

0^ G L2(0,T;L2(R+)), N{T)<m 

inf l(v^ + 0),(e + c)}(t,O>i«- 

[0,T]xR+ 



(3.6) 



(3.7) 



for some positive constants m, m. 

Proposition 3.1 (Local existence) Let (0o,V'o,Co) ^ i/"^(R^). If || (0o, ^Ao, Co) ||i < m and 
inf [o,T]xR+{(^ + 0), (0 + C)}(^, — Zli, then there exist 5i and to = tQ{rn,fn) > such 
that if the wave amplitude satisfies 6^ + 6'~^^ < 61, then the half space problem (13.51) 
admits a unique solution (0, C)(t, ^) G Xs,,2m{^yto)- 

To prove Theorem 2.1, it is sufficient to prove the following a priori estimate. 

Proposition 3.2 (A priori estimate) Suppose that the half space problem (13. 5p has a 
solution (0, ?/^, C)(t, ^) G X^^^^[0,T] for a suitably small constant £0 > 0. There exists a 
positive constant ^2, such that if the wave amplitude satisfies 6^ + 5*"^ < ^2, then the 
solution (0, ■?/',C)(^,0 satisfy that for V t G [0,T], 

ft 



N{t) + 



r, 



+ 



Q)(r,-)||?t/r<C(iV(0)+(52 + eB), 



(3.^ 



where the positive constant C is independent of t. 
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3.2 Boundary estimates 



In this section we will obtain the boundary estimates needed in the analysis below. From 
the definition of the viscous contact wave fl2.19p and its property f l2.2ip . we have the 
following estimates, which are very important in the boundary estimates, 



(t;, -V^^,M,-f/^^,^^, -0^^)(t,O = 0(1)5^^6-^"^^^^, as |e + (T_t|^oo. (3.9) 



So on the boundary ^ = 0, 



. V', C) {t, e = 0) =iv^- V, M_ -U,9^- 6) (t, e = 0) 

= {v, - y^^, - f/^^, 9, - e^^)(t, ^ = 0) 

= 0(l)(5^^e-"i*, as t ^ +oo. 



(3.10) 



Thus we have the following lemma. 

Lemma 3.1 (Boundary estimates) There exists the positive constant C such that for any 
t > 0, 



"*|(0,^,C)|V,e = O)rfr<C7(O', 



% ®5nCi/ c — r\\^^ ^ f ni/- l|2 I »/■ ||2\ 7 I ^ (xCD\2 



V V U Jq 
\-^i%' + ^^]{r,^ = 0)dT<u f\M'dT + C^ fu,fdr + C{5^^)\ 

V V Jq Jq 

(3.11) 

V 

where w = — , z/ is a positive small constant to be determined later and C^, is a positive 
constant depending on u. 

Proof: The proof of (13.1 II) i is a direct consequence of fl3.10p . 
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Now we prove flS.llI) ?. 

<c fmm + \UiM){r.i = mr 
t t 

<C [ \^^J^\mT,^ = 0)dr + C [ (|^p + |0|2)(r,e = O)cir 



<C f sup \^^{T,O\-mt,^ = 0)\dr + C{6'^y 
Jo 5e[o,+oo) 

<c f m^, oii^ ■ WMr, oii^ ■ mr,^ = o)\dT + c{6^^r 

Jo 

<v f\\\^^^[r,.)r+\\Urr)\?)dr + C. f \i,{T,i = Q)\'dr + C{5' 
Jo Jo 

< ^ ■)r+ mr, ■)r)dr + 



CD\2 



So the proof of flS.lip o is completed. Similarly, we can obtain (IS.lip cj. 
Then we will verify the inequality fl3.11l) ^. Notice that 

0^ V^(j) 



and 



So we have 



V V V V vV ' 

Vt Vt Vt a^v^ + a_V^ + 
V V _ V V V 

V V V 



*[-^^(|)= + 4](r,5 = 0M. 

t 

2 I l„/. |2 I I J|2 I L;,|2\ 



<C [ (|0^|2+|^^|2+|0|2 + |^|2)(^^^^o)rfr 

Jo 

<C f sup |^5(r,Opt^r + C(5^^)2 
Jo Ce[o,+oo) 

<C f\\Ur,-)\\m{r,-)\\dr + C{6^''? 
Jo 

<^ twAdrrWdr + C. T ||^^(r, Of c/r + 
Jo Jo 
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(3.12) 



where in the second inequahty we have used the fact 
Now Lemma 3.1 is proved. 



3.3 Energy estimates 



In this section we will prove the a priori estimate in Proposition 3.2. Firstly we have the 
following Lemma: 

Lemma 3.2 There exist a constant C > such that if the wave amplitudes 6^ , 5^^ and 
the constants e, Eq are small enough, then we have V t G [0, T], 

V^,C,0c)(t,-)ll'+ r||(05,^5,a)(r,-)ll'^^+ f [ U^{<l>' + e)d^dT 

Jo Jo JR+ 

< C||(0o,^o,Co,0o€)f + + + Al + r)-^||(0,^,C)(-,r)fdr + l 




JR+ 



Jo 

Proof: Let 

(^{z) = z - 1 - Inz. 
Similar in we can get the following estimate 



-Cn(? + CT_T)^ 



C 



, ^^^^ 



(3.13) 



(3.14) 



where 



i/i(«,0 = -<T_/,(«,{) + (p-PW'-A.( 



V V V V 



+ 



(3.15) 
(3.16) 



+$( 



9 



V V ^9 



v9 



v9^ 



V V 9^ 



i_ Kce| 

1/^ 9^ 



(3.17) 
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Note that 



$(1) = $'(1) = 0, <l>"{z) = z^^ > 0. 
So there exists a positive constant C such that 

C-\z-l)^ <^z) <C{z-iy, 



if z is near 1. 

Using the a priori assumptions N{T) < sq for suitably small constant Eq, we can get 



and 



c-vp < H^) < m', c-'\c\' < $(|) < cicr, 



c-'\i<p,C)\'<H^^)+iH^)<c\{<p,0\'. 



Substituting ( 13.180 and ( 13.190 into ( I3.17P and using Cauchy inequality imply 



Q3< 



+ 



Ave AvO'^ 

+(|(i^f,f/f,ef)|Meg|) + |g 



(3.18) 



(3.19) 



(3.20) 



By the fact 



we have 



1/(01 = 1/(0)+ / /.rfi/l< 1/(0)1 + v^||/dl, 



f [ \{vf, f/f , ef , eg) I (02 + e)d^dT 

Jo J R+ 



< C6^ 




'0 JR+ 



C)P(r,O)+e||(0^,a)f) did 



T 



<C6^ \{<P,Onr,0)dT + C6^ mXdrdr 
Jo , Jo 



Jo 

By the properties of the viscous contact wave, we can obtain 



(3.21) 



(3.22) 




'0 J R+ Jo Jlt+ 

(3.23) 
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Using the definition of tlie approximate rarefaction wave, we have 

r / (|(^f,f/f,ef)|Me||)(</>2 + m 

-'0 J R+ 

< r(ll(n^f/f,0f)f + l|0glUOII(0,C)lliocrfr 

<Ce'^ /"'(l + ^)-^ll(0,C)llll(0«,a)IMr 



(3.24) 



<Ce-^ / (l + r)-^||(0,C)frfr + C5B / \\{<t>^,Cd\?dT 
Jo Jo 

where in the second inequahty we have used 

||(ef,yf,f/f)f <Cei(l + t)-i, 

and 



||(eg,i^|,f/5)lki<C5t(i + t) 



13 
' 16 



if we let g > 16 in Lemma 2.3. 

Now we estimate the terms Qiip, on the right-hand side of fl3.14p and the term 
|Q2|(0^ + C^) on the right-hand side of fl3.20p . Due to the estimation of Qi in (13.31) . we 
have 

r / \Q,mdr <C fuWLAQiWLldr 
Jo Jn+ Jo, * 



< C 



dr 



3(1 + r) 12 (ir 



\l+T)-^dT+l 



(3.25) 

Similarly we can calculate the term and |Q2|(0^ + C^)- 

Integrating (I3.15P over R"*" x [0,t\ and using the boundary estimates in Lemma 3.1, 
we can obtain 



ll(0,^,C)(t,-)lP 



\m,Q){r,-)rdT+ 







f/f(0^ 



C')d^dT 



< C||(0o,V^o,Co)f + C(5^ + 5^^ + ^^) / (l + r)-Tt||(0,V',C)f^r + l 



+Cu / ||(V^55,Cd(r,-)rc?r + C<5' 



CD 




{l + T)-'e ^T^\{^,0\'d^dT. 



JR.^ 



(3.26) 
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Now we estimate H^^H^- Let 



v 



From the system (13.51) 9. we have 



Multiplying the above equation by ^ and noticing that 

Re RQ 



ip-Ph 



V V 



+ {p-p)^-Re^{ 



V v' 



we can get 



2 ^v^ 



V V 



RQ 



-ip-P)^ + RQ^i 



V V' 



Qi 



!1 

V 



V - ' V 

Integrating the above equality and using the boundary estimate (13.111) . we obtain 



< 



R+ 



R+ 



-c 



^(^)2-^^' 








RB_,v^^2 
2v ^ V 



d^dr 



(0,^)d^ + CiS^y + u I l|^55(r,-)frfr 



(3.27) 



dT + C 








Using the equality (I3.12p . we can get 

c-\m' - m') < Cfr < cm' + m')- 

By the estimation of Qi in (13. 3p . we have 

ri|gifc^r<C((5^ + (5^^ + £i). 
Jo 



(3.28) 



(3.29) 



Similar to (I3.22I) - (I3.24I) . we can compute the last term in the right hand side of (13.271) . 
Thus we can obtain 

ft 



'dr < C\\{^o,i^o,M)f + C\\{<P,m,-)\\ 



+Cu / U^^{T,-)rdT + C5' 



CD 




(l + r)-ie ^T^\{<PX)?didT 



(3.30) 



JR+ 



+C 



'5,C5)frfr + C(5^ + 5^^ + et) / (l + r)-S||(0,^,C)||^dr + l 
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Multiplying the inequality fl3.26p by a large constant Ci > 0, and adding it to fl3.30p . we 
complete the proof of Lemma 3.2. 

Now we derive the higher order estimates. Multiplying the equation (13.5^ 9 by — 
we can get 



f)t-[M^-^h+^^^ = ^^^v^^^^+{ip-p)^-^^[u^i-^ (3.31) 

From the boundary estimate 







t 

2 I l„/. n\\2\ 



< C / i\MT,0)\' + \Mr,0)\')dT 







< c m\MdT+c{6' 



CD\2 



< V 



Jo Jo 




we can get the following inequality by integrating (13.311) over R+ x (0, t) 

IM'it) + f \m?dr < CU,^r + C{6^^f + C.^ Al + r)-i ||(</., Ofdr 
Jo ^ Jo 

+C II (05, Cc) frfr + / / {l + r)-'e ^^UX)?didr, 

Jo Jo Jii+ 

(3.32) 

where we have used the following estimation 

i^^ii^^ii^^^iderfr <c fmrndmUfdr 

'0 ^0^ 

<c f ll^dlll^^dl^ll^dl^^^ 

J\ t 

<u [ ii^5dr^^+c^supii0di' / ii^di'^^ 

t * t 

<u [ \M'dT + C,e', [ Wi^^fdr. 
Jo Jo 

Multiplying ([33])3 by — C^c, almost similar to the estimates for ||'?/'5p(t), we can obtain 

'it)+ r||C«frfr<C||CodP + C(0' + C^^^ [\l + r)-h\{<l>X)rdT 
Jo ^ Jo 

+C / II (05,^5, C5)frfr + C(0' / / {l + r)-'e ^^|(0, Oprf^rfr. 

Jo Jo Jn.+ 

(3.33) 
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Combining Lemma 3.2 and the higher order estimations (13.321) and fl3.33p . we have the 
following Lemma: 

Lemma 3.3 If the wave amplitudes 5^ , 5'"^ and the constants e, Sq are small enough, 
then we have V t G [0, T], 

ii(0,^,c)(t,-)ii?+ rii0di' + 11(^^^)11?^^+ f I uf{<p' + e)didT 

Jo ^ Jo JR+ 

< C||(0o,V^o,Co)||? + C^(5^ + 5^^ + ^^)f / (l + r)-i||(0,^,C)irrfr+ll (3.34) 



+C5' 



CD 



-, -co(e + cT_T)^ 



JR 



In order to close the estimate, we only need to control the last term in (I3.34p . which 
comes from the viscous contact wave. So we will use the estimation on the heat kernel in 
[5] and [7]. 

Lemma 3.4 Suppose that h{t,C,) satisfies 

he L°^{0,T-L\R+)), h^e L\0,T;L\R+)), ht - (x_h^ e L\0,T; H-\R+)), 
Then 



1 + ry^h-^e ^ d^dr 



JR+ 



where 



IIM0,Oir+ / h\r,i = Q)dr+ / \\h\\''dr+ / {ht - a_h^,hgl) u^^H-^dr 



(3.35) 



g0{t,i) = -{l + t)- 



e i+* di], 



and /3 > is the constant to be determined. 

The proof of Lemma 3.4 can be done similarly in [3J. The only difference is that the 
space we considered here is on the half line and the boundary terms should be treated. 

Lemma 3.5 There exist a constant C > such that if 5^^ and Sq are small enough, then 
we have 



co(g + CT_r)-' 

1 + T 







< c 



1 + r 



-\{(t^,^X)?didT 



(5« + sc^ + e-^) + II (00, ^0, Co) ir + II (0, ^, C)(t, ■) \[ 



(3.36) 
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Proof: From the equation fl3.5p 9 and the fact p — P = BL^i.^ have 



V V V 



Then we get 



{RC - P0)5 = ^^^-^{V^ + 0^) - v{i^t - + livC^ - - vQ,. (3.37) 



Let 



G,(t,0 = -(l + t)"' / e~dri 



+00 



where a is a positive constant to be determined later. Multiplying the equation (I3.37P by 
G„(i?C - P0) gives 



G^{RC-Pcl>f] , [EC - Pct>f 

- ^ 



G^{RC-P(t>f 



(3.38) 



+/iG'„t;(i?C - P^)(!^ - ^)5 - Go.v{Ri - P<P)Qi. 
V V 

Note that 

-Gc.v{RC - P0)(^i - a_z^5) = -[G^v{RC - Pct>)i,]t + [G^v{RC - Pct>)^]^ 



and 



+ [(G,t;)t - a_{G^v)^]{RC - P0)^ + G^vi,[{RC - P^)t - a^RC - P^)^], 

(PC-P+0)*-^-(i?C-P0)c 

= (PC - Rs-Q) - {Pt - ^-Pe)^ - P(0i - ^-0?) 

t;\/t>l/ . 



(3.39) 



-(P, - a_P5)0. 
And using the equality 



- G^v^Pi/j^^ = -[^G^vP^]^ + ^vPiG^)^^ + i{vP)^^, 



(3.40) 
(3.41) 



we can get 

a{« + <T_t)2 



e 1+* 
2(1 + t) 



■[(PC - P0)' + iPv^p^] = [Go.v{RC - P(P)tP]t + H2dt, + Qa, (3.42) 
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where 



V V 



V V 



and 



Q4 



-[{G^)t - a4G^)^]v{RC - P0)^ - G^u^iRC - P0)^ 
+(7 - l)GaV^ 



Note that 

||G,(t,-)l|L- <C„(l + t)- 
thus integrating (13.421) over R+ x (0, t) gives 




JR^ 



1 + r 



.[{RC-Pct)f + i,^]didr 



< C 



C / (1 + r) 



^,C)(r,-)fc^r + C / ||(05,^5,a)(r,-)frfr 



/ ||(^^^,C5^)(r,-)rcir + CT 



CD 








n(S + g--T)^ 

e 1+^ 
1 + r 



|(0,C)PrfCc^r. 



(3.43) 



(3.44) 



+/.[G„t;(PC - Pm"^ - ^) + (7 - l)«:(G„t;^)«(^ - ^) 
+GaviRC - P(p)Qi + Gc,v^{Pt - a.P^)<P, 

From the boundary estimate (13. lip , we have 



iJ2(r,e = 0)rfr<C(O' + ^ / ||(^«,C?,^cc,C«)(r,-)rrfr. (3.45) 



(3.46) 



In order to get the desired estimate in Lemma 3.5, we must derive the other similar 
estimates from the energy equation ([33])3. Set 

= PC + (7 - l)P+0 

in Lemma 3.4. Thus we only need to compute the last term in (I3.35p . From the energy 
equation ([3]5])3, we have 

h, - a.h^ = {Pt - a.P^)<P - (p - P)u^ + «:(^ - -^)^ + " if) " Q2, (3.47) 

V V V V 
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Thus 



/ / [[P,-a^P^)ct>-{p-P)U^]hgldidT+ f I {p-P)^^hgldCdi 
Jo , Jo Jn.+ 



+ / ['^ 

'0 



+ 




4_0i 

■V V 
^1 



6 



Jr+ V 



)hgl]{T,^ = 0)d 



T+ I I - ^){hgp)^d^dT 





JR+ V V 
Q2hgld^dT 



(3.48) 







i=l 

Note that 

\\g(3{t,-)\\L^<C(,, 

we can estimate Ji{i = 1,3,4,5,6) directly. In order to estimate J2, from the mass 
equation fl3.5P i . we have 

{p-P)^^hgl 
= hgp{(f)t - a^(f)^) 



V " 2v 

2{j-l)^h'g}-jPh4>'gp 



2v 

-fPh<f - 2(7 - l)h^<i) 



,2^2 



2{^-l)<Ph^gl~^Ph<P^gl 



+ 



gi3[{g(3)t-(y-{9i3)(\ 



2v 

-fPh<i? - 2(7 - 1)/IV 



2(7 - l)glct>h , ^Pglct>^ 



+ 



2v 



{ht - a^h^) + 



2v 



Now each term can be estimated directly, the detailed proof can be seen in [3]. Remark 
that here we need to compute the boundary terms. Therefore, taking the constant /3 = y, 
we can get from Lemma 3.4 that 



co(e+CT_T)^ 



'0 + ^ 




-hHidr < C (5^ + 5^^ + el) + l|(0o, V^o, Co)!!' + !! (0, V^, C)(t, ■) !!' 




CD 











+C(<5^^ + £o) 




CO{e + <T_T)^ 

{l + r)-'e ^T^\{^,0\'d^dT. 



JR+ 



(3.49) 
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Taking a = cq in fl3.46p and combining the estimates fl3.46p with fl3.49p yield the desired 
estimation in Lemma 3.5 if we choose suitably small constants 5'"'^ and Eq. 

Now from Lemma 3.3 and Lemma 3.5, if the wave amplitude 6'"^ and the constant i' 
are suitably small, we can get 



t 

2 



ii(0,^,c)(t,-)iif + m' + \mx^mdT<c\\{<j>o,^oXo)\\i 







(l + r)-i5||(0,^,C)|pt/r + l 







Finally, Gronwall inequality gives the a priori estimate in Proposition 3.2: 







Thus we complete the proof of Theorem 2.2. 

The proof of Theorem 2.1 can be done along the same line as Theorem 2.2, we omit 
it for brevity. 

Appendix: Proof of Lemma 2.1 

Now we give the rigorous proof of Lemma 2.1. Firstly from (12. 9p . u+ > 0. Thus if < 0, 
then there is no solution to (12.101) or (12.131) . Now assume that > 0. Then we can 
compute that the determinant of the matrix J defined in (12.120 

K/i(7 - 1) fi;/i(7 - 1) 

So we can divide it into three cases according to the sign of the quantity — 1. 
Case I (Supersonic): M+ > 1, then det J > 0. We can easily know that J has two 
positive eigenvalues. Thus the ODE system (I2.13P has no solution. 

Case II (Transonic): M+ = 1, then det J = 0. One of the eigenvalues of the matrix 
J is zero, the other one is positive. This case is a little subtle. Firstly we can choose a 
nonsingular matrix P such that P~^JP changes into a standard form. For example, let 



P 



( k(7-1)^ fi«^7(7-l) 
i?K/i+K(7— 1)^ [i?K/^+K(7— 



M"4 



«;(7-l) 



then 



P-'JP Q 1 := A, 
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where Xj is the positive eigenvalue of J given by 



7 



Let 



W 



Wi 

W2 



+ 



R 



«:(7 - 1) 



■U+ > 0. 



P 



QB 



we have 



where 



and 



We can rewrite fl2.13p as 



W^=AjW + G{W), 

G{W) = P-^F{PW) 
F(PW)- ( ^'^^^^ 



W^^ = XjWi + Gi{Wi,W2] 

W2^ = G2{W,,W2), 



(A.2) 



(A.3) 



(A.4) 



Obviously, there exists a suitably small neighborhood ^^^^(0, 0) such that {Gi, G2)iW) 
is analytic. And in this neighborhood, if \Wi\ <C \W2\, then 



G2iW) 



i?VK(7-l)(7 + l) 
■2[i?7/i + /€(7-i; 



212 



+ o{Wi 



(A.5) 



From the geometric theory of the automatons ordinary differential systems, we know that 
the equilibrium state (0,0) is a saddle- node point to the system flA.4l) . And (0,0) is an 
attractor whose trajectory, denoted by F, is unique and tangent to iy2-axis at (0,0). From 
the uniqueness of the attractor trajectory F, we know that only when (VFi,VF2)(0) G F, 
there exists a solution to (lA.4p . otherwise, there does not exist solution to fl2.13p . When 
(W^i,W^2)(0) e F, the solution {Wi,W2){i) satisfy that \Wi{i)\ < 1^2(01 if ^ is large 
enough, thus we have 



(A.6) 



where < ai < o"2 are two constants. 
So we can get 



\{w^,w2m\ < c 



1 + 6^^ 



(A.7) 
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where = \{Wi,W2){0)\ = 0(1)|(m+ - 0+ - ^_)| is small enough. 

From (]A.2p . we can get the BL-solution {U^ ,<d^){C,) in the transonic case (M+ = 1^ 
satisfy that 



fJ,U+ 


k(7 - 




POO 


u 




h 





tB 



u4U^r , (2-7)«H 



k(7 - 1) 



27K 



B\2 



+ 



^ UBqB_J_^^B^3 

k(7-1) 2k^ ' 



d^, (A.8) 



Case III (Subsonic): < 1, then det J < 0. One can see that J has one positive 
and one negative eigenvalues. Similar to Case II, we can choose a nonsingular matrix P 
such that P-^JP is in a standard form. Now we give the detailed procedure for choosing 
the matrix P. Firstly, let 

1 

ai 1 

where the constant ai is to be determined, then 



Pi 



Pi 



-1 



Assume that 



Pi'JPi 



1 
-ai 1 



mil nii2 
m22 



an, 



(A.9) 



where the constants rrin, 17112, will be fixed when Oi is determined. From m2i = in 
( lA.gp . we get a equation of ai: 



R 



al + 



I.e., 



ai 



+ 



Ml-i - 1 



Ru. 



«:(7-l) 



ai — 



/i 



fli 



0, 



M2i?7 k(7 -l)J u+ MlR-^n 
Then we can solve the above equation to obtain 



Oi = CiU+ < or ai = C2M+ > 0, 



where Ci < min |0, 
the following equation 



, C2 > max 



k(7-1)' k(7-1) MlR-r 



M|7- 1 
MlR-f ~ k(7 - 1) 



MiR-fK 



(A.IO) 

> are the solutions of 
(A.ll) 
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Without loss of generality, we choose ai = 02^+, then we can compute that the matrix 
on in ([Oil 



A^, 



where Xj = ^ ^^2^^ + £2^^+ > and Aj = i^-^r^^zx) ~ ^•1^'^+ < are the two eigenvalues 
of the matrix J. 



Then we can choose a matrix 



such that 

Then we can get 

Now we set 

Then 
Let 

we have 

where 
and 



1 a2 
1 



^2 



P^-igjlP^ = ( ^2 ) := Aj, (A.12) 



^2 = JT\ T^- (A.13) 



\ 21. 

u+ 

C2U+ 1 + a2C2 



P = P1P2 
p-'jP = Aj = dmg{\),\]}. 



W={ wl]--=P-'[ Is ), (A.14) 



W): = KjW + G{W), (A.15) 
G{W) = P-^F{PW), 



We can rewrite (1A.15|) as 



= X]Wi + Gi{Wi,W2) 

= X]W2 + G2{Wi,W2) 



(A.16) 
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From above, one can easily know that Gi,G2 are analytic with respect to (1^1,1^2) 
near (0, 0), then the equilibrium point (0, 0) is the saddle point of (]A.16p . i.e., in a suitably 
small neighborhood f^^(,(0, 0), there exist two opposite attractor trajectories Fi, F2 tangent 
to PF2-axis at (0,0). Let = Fi U F2, then is a center-stable manifold. Only when 
(H^i,W2)(0) G tM, there exists a solution the ODE system ( 1A.16I1 . In such case, there 
exist two positive constants cr3,o"4 which is close to — such that 

- a3W2 <W2< -(TaW2, i E R+. (A.17) 

So we have that there exist positive constants c and C such that 

\iW^,W2m\<C5^e-'^, eeM+, (A.18) 

where 5^ = \iWi,W2)iO)\ = 0(1)|(m+ ^+ - is the amplitude of the BL-solution. 
The BL-solution ([7^,6^) satisfies 



tB 



(1 + a2C2U+){U 

1 + a2C2 



U- 



02(6^ - e 



B\2 



a2 
u 



+ 



+ 



R 



K{-i - 1) 



2k 



B\3 



(A.19) 



Now we complete the proof of Lemma 2.1. 
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